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The Bogoliubov - de Gennes equations are solved for the Coulomb Bose gas describing a fluid 
of charged bosons at finite temperature. The approach is applicable in the weak coupling regime 
and the extent of its quantitative usefulness is tested in the three-dimensional fluid, for which 
diffusion Monte Carlo data are available on the condensate fraction at zero temperature. The one- 
body density matrix is then evaluated by the same approach for the two-dimensional fluid with 
e 2 /r interactions, to demonstrate the presence of a quasi-condensate from its power-law decay with 
increasing distance and to evaluate the superfluid fraction as a function of temperature at weak 
■ coupling. 



o 



(N 
O 

Q 



PACS numbers: 05.30.Jp, Boson systems 



I. INTRODUCTION 



The fluid of pointlike spinless charged bosons embedded in a uniform neutralizing background has attracted attention 
in the literature mainly as a model in quantum statistical mechanics, which is complementary to the physically more 
relevant fermionic gas of electronstH. It was proposed by Schafrotbo as a model for superconductors prior to the BCS 
£ — ' theory and has received renewed interest after the discovery of ceramic superconductors! In some viewpoints a Bosc- 
Einstein condensate of .tightly bound pairs of small polarons could be a relevant model for high-T c superconductivity 
■ in the layered cupiatesu. The model also has some astrophysical relevance in describing pressure-ionized helium in 
cold stellar matteiulj and the fusion of a-particles inside a dense helium plasmaou. 

A number of theoretical and computational studies have been addressed to the three-dimensional charged-boson fluid 
(3D-CBF) at zero temperature. The properties of main interest for the theory have been the ground-state energy and 
structure and the static and dynamic dielectric response. The early theoretical work was concerafid-adth evaluating the 

aground-state energy and the elementary excitations in Ihc-weak coupling (high density) Iimitt3il3. Both variational 
calculations based on Jastrow-Feenberg wave functionsEjll3 and self-consistent treatments of correlationstjci have 
subsequently teen used to evaluate the intermediate and strong coupling regime. Quantal Monte Carlo studies of the 
3D-CBFE3l3l3 have covered the whole range of coupling strength up to the regime of Wigner crystallization driven 
by the Coulomb repulsions. Extensive data on the condensate fraction and the momentum distribution in dependence 
of the coupling strength have become available through the diffusion Monte Carlo (DMC) work of Moroni et al. 2 A. 

The properties of the 2D-CBF at zer-CL-temperature have also been investigated with both e 2 /r and ln(r) interactions. 
In the latter case Magro and Ceperleyo have shown, using a sum rule argument from Pitaevskii and StringariEII, that 
the presence of the long-wavelength plasmon in the excitation spectrum rules out the existence of a condensate even 
at zero temperature. They further showed from a DMC study that the noncondensed fluid exhibits a power-law 
decay of the one-body density matrix and that above a threshold density the momentum distribution diverges at low 
momenta, although no condensate forms. 

In the 2D-CBF with e 2 /r interactions (henceforth referred to as quasi-2D or Q2D-CBF) a condensate is present in 
the, ground state, but its density vanishes at finite temperature. This behaviour parallels that of the neutral 2D Bose 
gastEro. The theory p£ correlations in the latter system has been developed by Kagan, Svistunov and ShlyapnikovB3 
(see also Kagan et alSS). As its temperature is lowered at constant density (or as its density is increased at constant 
temperature) across the Kosterlitz-Thouless transition, a weakly interacting gas enters the superfluid regime in which 
the phase correlation length R c is much larger than the density correlation length r c . In this situation the one-body 
density matrix p(r) decays asymptotically to zero with an inverse-power law rather than exponentiallyoa. The 
idea of a "quasicondensate" emerges from the behaviour of p{r) at intermediate distances r c <C r <C R c . The local 
properties of the quasicondensate are the same as those of a genuine condensate, so that it turns to the latter as the 
phase correlation length starts to exceed the dimensions of the sample. 

The main purpose of the present work is to study this behaviour in the Q2D-CBF, using the Bogoliubov approach 
to describe the charged fluid both in the ground state and at finite temperature in the weak coupling regime (corre- 
sponding in this case to high density). The Bogoliubov - de Gennes equations are presented for convenience in Sect. 2 
and are first solved for the 3D-CBF in Sect. 3, where the approach is quantitatively tested at T = by comparing 
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its results for the condensate fraction with the available DMC data. Section 4 reports our main results regarding 
the quasicondensate and the superfluid fraction in the Q2D-CBF. Finally, Sect. 5 gives a brief summary and our 
conclusions. 



II. BOGOLIUBOV APPROACH TO A CHARGED-BOSON FLUID 



The fluid of charged bosons on a uniform neutralizing background is described by the Hamiltonian 
H = J dr^{v) -I*) il>{*) + \Jdv J dv' ^ {v)^{t')V{\v - r'|)t/;(r')^(r) , 



(1) 



where V(r) — e 2 /r, ip(r) is the field operator and \x the chemical potential. The role of the background is to set to 
zero the average potential felt by each particle. The coupling strength is measured by the dimensionless parameter 
r s , defined by r s as — (4-7m/3) -1 / 3 in 3D and by r s aB — (nn)^ 1 / 2 in 2D with as the Bohr radius and n the mean 
particle density. _ 

We consider first the fluid at T = 0. The Bogoliubov transformation^ (for recent presentations seoo) introduces 
a macroscopic order parameter ip$ by writing the field operator as ^(r) = ipo + VK 1 ")- The operator ip(r) describes 
the gas of Bose particles promoted out of the condensate. This gas is treated in the Hartree approximation, assuming 
that the contribution from terms nonlinear in ip(r) is small. 

One finds fi = and, with the linear transformation 



V3(r, t) = Mr, t)a„ + <(r, t)4] 



(2) 



for the Heisenberg field operator, one has to solve the coupled linear equations 



.du„(r, t) 

'' dt 



1 

2m 



V 2 u u (r,t) + n / dr'V(\r-r'\)[u„(r,t)+u u (r',t) + v l/ {r',t)] 



(3) 



and 



. dv u (r,t) 
1 dt 



2m 



V 2 v„{r,t)+n / dr'V{\r-r'\)[v u {r,t)+v u (r',t)+u v (r',t)] 



(4) 



Here, no = V'o i s the uniform condensate density. The subsidiary condition 



]T Mr, *)<(r', t) - v v {v', t)v* u (r, t)] = S(r - r') 



(5) 



embodies the commutation rules on the field operators. 

In a uniform fluid the state index v is the wave vector k. Equations (3) - (5) are solved by taking Uk(r, t) 
Mk exp[i(k • r — e^t)] and Vk(r, t) = i>k exp[i(k ■ r — £ki)], with the results 



n V(k) 



,2 1 



2m 



(6) 



and 



vl = --{1-eZ 1 



n V{k) 



2m 



(7) 



where 



The condensate density is given by 



£k = 



n k 2 V(k) 


(-Y 
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\2m J 



1/2 



no = n — (^ t (r)i/'(r)} = n — v 



(8) 



(9) 
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and the one-body density matrix is 

p(r) = n + {ft (r)4>(0)) = n + ^ w£ exp(ik • r) . (10) 

We have assumed unitary volume in writing Eqs. (9) and (10). 

After these transformations and approximations the Hamiltonian (1) has been reduced to that of a system of inde- 
pendent bosonic excitations described by the operators au and a^. The extension of the theory to finite temperature 
is then effected by means of the Bose-Einstein distribution function 

h = [exp(/3e k ) - If 1 . (11) 

In particular, Eqs (9) and (10) become 

n =n-J2K + fk« + vl)] (12) 

and 

p(r) = n + [«k + M u l + V D] ex P(* k ' r ) » ( 13 ) 

k#0 

respectively. 



III. THE THREE-DIMENSIONAL CHARGED-BOSON FLUID 



We introduce the parameter A — (3no/nr^) 1 / 2 e 2 /(aBkBT) and the reduced distance R = (12no/nrg) 1 / 4 (r/aB 



3U/4/ 



Then Eqs. (9), (12) and (13) can be rewritten as follows: 



1 _ (H° 



-3/4' 



T=0 



2 1 / 2 r s 3/4 

3V*7T 



dx 



g( x ) 



2x 2 



(14) 



V re / In 



and 



p(r) 



2 1 / 2 / 
3!/% I 



-3/4 



n r s 



2 i/2 r 3/4 r dx if (l] 



T^O 



3/4 



3!/% 



exp[Ag(^)] - 1 



2a;' 



sin(i?a;) 
Rx 



i + 



exp[ J 4g(a;)] — 1 



2x' 



(15) 



(16) 



where we have defined f(x) = 1 + 2x 4 and g(x) = (l + x 4 ) 1 ' 2 . Of course, Eq. (14) can also be obtained from Eq. (15) 
in the limit T — > 0. In the following we use the energy e 2 /as as the unit of the thermal energy ksT. 

A numerical solution of Eqs. (14) and (15) for no/n in the physical range < no/n < 1 can be found for all values 
of the system parameters. We should bear in mind, however, that we are using a weak-coupling theory so that the 
results can be significant only when the depletion of the condensate is small (i.e. for no/n close to unity, in line of 
principle) . This statement is quantitatively tested in Table jjj, where we compare our results for no/n at T — with 
the DMC data of Moroni et alEB over the whole fluid range up to Wigner crystallization. It is clear from Table | that, 
rather surprisingly, the Bogoliubov approach is almost fully quantitative up to r s ~ 5, i.e. for values of no/n down 
to almost 0.5. Very similar results are obtained in the same range of r s by the integro-differential equations approach 
of Cherny and ShanenkcO. 

Table O reports our results for no/n as a function of temperature for r s = 1 and r s = 2. In the lack of data 
for a quantitative test, one may hope from the test shown in Table | that the Bogoliubov approach could again be 
reasonably accurate for values of no/n larger than 0.5. According to this crude criterion, it appears from Table [] 
that at such weak couplings the theory could perhaps be useful up to a fairly sizable value of the reduced temperature 
fcsTas/e 2 - perhaps as large as unity for r s — 1. 

Finally, the one-body density matrix p(r) (in units of the particle density n) is shown in Figure 1 for r s — 1 (left 
panel) and r s = 2 (right panel), at various values of the reduced temperature. The asymptotic value of p[r)/n in the 
limit r — > oo is the condensate fraction no/n, that we have already presented in Table ||. 
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IV. THE QUASI-TWO-DIMENSIONAL CHARGED-BOSON FLUID 



Taking V(k) = 2ire 2 /k, Eq. (9) for the condensate fraction at T = yields 



1 



no\ /no\ 
n I \ n ) 



-2/3* 



r 2/3 



T=0 



dx 



2x* 



(17) 



where we have defined f(x) = 1 + 2x 6 and g(x) — (1 + x 6 ) 1 / 2 . The numerical solution of Eq. (17) yields the values 
of the condensate fraction which are reported in Table III over a range of values for the coupling strength well below 
the Wigner phase transitioned at r s ~ 35. We should expect that the role of correlations becomes more important in 
lowered dimensionality and indeed the crude criterion introduced in Sect. 3 suggests that in the Q2D fluid at T = 
the Bogoliubov approach may be useful only up to r s ~ 1. 

It is easily seen that the condensate fraction must instead vanish at T / 0. Indeed, if we assume hq/ii / then 
Eq. (12) yields 



no\ /«£\~ 2 / 3 
n J \ n J 



,2/3 



dx 



i + 



cxp[ylxg(x)] — 1 



2x ; 



(18) 



where we have set A = 2{no/nr 2 ) 2 ' 3 e 2 /(asksT). The second term in the square bracket on the RHS of Eq. (18) 
yields a contribution of order a; -1 to the integrand for x — > 0, so that the integral diverges. On the other hand, the 
solution no/n = at T ^ is consistent with the Bogoliubov - de Gennes equations: in this case they yield = 
and Uk = 1, and the particle density is related to the (now finite) chemical potential by 



$Nexp 







2m 



- 1 



(19) 



We can at this point introduce the quasicondensate for the Q2D-CBF at finite temperature. The essential point is 
the powetrlaw decay of the one-body density matrix, which becomes slow at sufficiently low temperature and weak 
couplingE2rE3. We first show this by an heuristic procedure that we shall justify later below. We isolate the singular 
term in p(r) and resum it to infinite order to obtain 



p(r) = p{r) exp[-A(r)] , 



(20) 



where with the help of Eq. (13) we have set 

d 2 k 



A(r) 



2 [l-cos(k-r)]^/( £ - k ) 



= 2 



(27T) 

2 \ 1/3 roc 

h 



dx 



1 - J (x 2 R) 



g{x) {exp[Axg{x)} - 1 j 



(21) 



and 



m = 1 _ i 



d 2 k 



n J (2tt) 2 

> 2/3 

n r s " 



[1 - cos(k • r)] 



1 / n I/(fc) + fc 2 /2m _ A fc 2 /2m 



dx [1 - J (x 2 R)} 



: /(e k ) 



Ax b 



9( x ) g(x) (exp[Axg(x)] - l) 



-2x 6 



the quasicondensate density ho being defined by 



no = km p(r) . 



(22) 



(23) 



In these equations e k = {h V(k)k 2 /m + k 4 /Am 2 ) 1 ^ 2 , A = 2(n /w 2 ) 2 / 3 e 2 /(a B k B T), R = 2(h /nr 2 ) 1 / 3 (r/a B ) and 
Jo(y) is the Bessel function of zero-th order. 

Figure 2 (upper panels) reports our numerical results for p(r) at two values of the coupling strength and at various 
values of the reduced temperature. It is evident that at low temperature the decay of the density matrix in space is 
slow and the notion of a quasicondensate thereby acquires physical significance. 
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A power-law decay of the one-body density matrix at low temperature and coupling strength can be demonstrated 
directly from Eqs. (20) - (23). The function Jq(x 2 R) in the integrand in Eq. (21) provides a lower limit of integration 
going as r -1 / 2 , wiule the upper limit is set by a cut-off wavevector fc = 1/L associated with the quasicondensate 
region (see PopovEj). Eq. (21) thus yields 



-1/2 



A(r) -> 2i- 1 (w 2 /n ) 1/3 / dx\ = 2A-\nr 2 Jh ) 1 ^ \n{r/L) 1 '' 1 . (24) 
Hence, since p(r) — ► no we find from Eq. (20) 



-1/2 



with the value of the exponent given by 



p(r) -» Mr/L)~ a (25) 



-=^ s ^-. (26) 



Figure || (bottom panels) evidentiates the power-law decay of the density matrix by plotting the curves in log-log 
scale, as compared with the aympotic behaviour predicted by Eqs. (25) and (26) (dots in the Figure). 

A. Asymptotic behaviour of the single-particle Green's function 

We can actually show that the power-law decay of p{r) that we have obtained just above derives from the correlations 
between phase fluctuations in the Q2D=CBF. We follow the method proposed for the neutral 2D gas in the work of 
Popo\c3 (see also Fisher and Hohenbergo) . The fluctuations of the phase 0(x, t) determine the single-particle Green's 
function in the low-momentum regime (below the cut-off momentum fco) according to 

G(x,t;xi,ti) ~n exp j-~([^(x,r) - 0(xi, n)] 2 ) 1 . (27) 



From Eq. (19.16) in Chapter 6 of Popov's bookO we find 

V{k) 
uj 2 + fioV(k)k 2 /m 

for the Q2D-CBF at long wavelengths and frequencies, so that 



<*(k, uM-k, -»)) - .. 2M ">, (28) 



V(k) 



2(3 ^ lu 2 + fi V(k)k 2 /m 

' k<k u v ' ' 

x | exp[i(k • x — lot)] — exp[i(k • Xi — ujt{)\ \ 2 . (29) 

In Eq. (29) we carry out the summation over the Matsubara frequencies for t\ = r + and keep the most diverging 
contribution in the integral over the momenta to obtain 

i([0(x,r) -0(x 1 ,r+)] 2 ) ^aln^+const. (30) 

where r — |x — xi| and the quantity a is given by Eq. (26). Using this result in Eq. (27) we see that the Green's 
function decays to zero with the law r~ Q . 

B. Superfiuid fraction and quasicondensate fraction 

As already noted in Sect. 1, the notion of a quasicondensate becomes meaningful at temperatures below the 
Kosterlitz-Ihouless transition. A superfiuid component is therefore present at these temperatures and its density n s 
is given bytj 

^ = l-^Vfc 2 6 ^ £k) 12 (31) 
n 2nm [exp(/3e k ) - l] 2 
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From Eq. (31) we have 



n 
n 



r s n 
n 



2/3 



dx ■ 



exp 



|exp A(x 



2-4)1/2 



i 2 ' 



(32) 



- 1 



The value of the superfLuid fraction n s /n that we obtain from Eq. (32) are reported in Table |v|for two values of the 
coupling strength and at various values of the reduced temperature. 

In Table IV we also show the values taken by the quasicondensate fraction no/n at the same values of the coupling 
strength. The quasicondensate density is similar to the superfiuid density at very low coupling, but rapidly decreases 
as the coupling strength is increased. 



V. SUMMARY AND OUTLOOK 



Summarizing, we have studied quasicondensation and superfluidity in a weakly interacting 2D fluid of charged bosons 
with-e 2 /r interactions at finite temperature. By comparison of results on the 3D fluid with diffusion Monte Carlo 
dataEj we have found that the Bogoliubov approach may yield quantitatively useful predictions over a surprisingly 
wide range of coupling strengths and of deviations of the condensate fraction from unity. At finite temperature the 
behaviour of the charged 2D gas is qualitatively wholly similar to that of its better known neutral analogue: well 
below the Kosterlitz-Thouless transition a slow power-law decay is seen in the one-body density matrix, heralding 
extended-range correlations in phase fluctuations and the formation of a condensate over regions of finite size. 

From the diffusion Monte Carlo results of Magro and CeperleyE^I it appears that the charged 2D fluid with loga- 
rithmic interactions may also show quasicondensation, even in the absence of macroscopic condensation in the ground 
state. Other interesting questions in this area regard how the transition from 3D to 2D behaviour is effected and how 
one may develop a sound theoretical description of the momentum distribution with increasing coupling strength. 
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TABLE I. Condensate fraction in the 3D-CBF at zpfio temperature from the Bogoliubov approach (B), compared with the 
diffusion Monte Carlo data (DMC, from Moroni et al£B). 



1 2 5 10 20 50 100 160 



(n /n) B 


0.818 


0.722 


0.549 


0.401 


0.264 


0.132 


0.072 


0.047 


(n /n) D MC 


0.827 


0.722 


0.542 


0.359 


0.206 


0.053 


0.0104 


0.004 



TABLE II. Condensate fraction in the 3D-CBF as a function of the reduced temperature from the Bogoliubov approach, for 
two values of the coupling strength r s . 



T 


{n /n)r s=1 


T 


(n /n) rs= 2 





0.818 





0.722 


0.3 


0.816 


0.1 


0.718 


0.5 


0.793 


0.2 


0.652 


0.75 


0.720 


0.3 


0.466 


1.0 


0.593 


0.32 


0.410 


1.25 


0.407 


0.34 


0.345 


1.4 


0.247 


0.36 


0.266 


1.482 


0.092 


0.38 


0.153 



TABLE III. Condensate fraction in the Q2D-CBF at zero temperature from the Bogoliubov approach. 



r a 


0.01 


0.1 


0.2 


0.4 


0.6 


0.8 


1 


2 


5 


n /n 


0.968 


0.863 


0.794 


0.700 


0.633 


0.580 


0.537 


0.398 


0.230 
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TABLE IV. Superfluid and quasicondensate fraction in the Q2D-CBF as a function of reduced temperature from the Bo- 
goliubov approach, for two values of the coupling strength r s . 



T (n 3 /w) rs = .i (n /n) rs=(hl T {n s /n) rs =i (h /n) rs = 1 

1 0.863 1 0.537 

10 0.991 0.861 0.2 0.995 0.536 

15 0.972 0.854 0.3 0.979 0.533 

20 0.944 0.843 0.5 0.905 0.517 

30 0.865 0.809 0.7 0.779 0.482 

40 0.763 0.758 0.8 0.696 0.457 




o 1 o 

02468 10 02468 10 

r/ a B r/ a B 

FIG. 1. The one-body density matrix p(r)/n as a function of r/as in the 3D-CBF. Left: for r s =l at values of the reduced 
temperature T equal to 0, 0.5, 1.0, 1.25, and 1.482 (from top to bottom). Right: for r s = 2 at T =0, 0.2, 0.3, 0.34 and 0.38. 
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FIG. 2. The one-body density matrix p(r)/n as a function of r/as in the Q2D-CBF, both in linear scale (upper panels) and 
in logarithmic scale (lower panels). Left panels: for r s =0.1 at values of the reduced temperature T equal to 0, 10, 15, 20, 30 
and 40 (from top to bottom). Right panels : for r s = 1 at T = 0, 0.2, 0.3, 0.5, 0.7 and 0.8. The dots in the lower panels 
give the asympotic behaviour of p(r)/n at large r as evaluated analytically in Eqs. (25) and (26). The length L in Eq. (25), 
which is not determined by the asymptotic calculation, appears in logarithmic scale as an additive constant and is here fixed 
by requiring that the analytic result overlaps the numerical one as r increases. 
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